Charge transfer via a two-strand superexchange bridge in DNA 
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Charge transfer in a DNA duplex chain is studied by constructing a system with virtual electrodes 
connected at the ends of each DNA strand. The system is described by the tight-binding model 
and its transport is analyzed by the transfer matrix method. The very weak distance dependence 
in a long (G:C)(T:A)>i (G:C)3 DNA chain observed in experiment [B. Giese, et ai, Nature 412, 
318 (2001)] is explained by a unistep two-strand superexchange bridge without the need for the 
multi-step thermally-induced hopping mechanism or the dephasing effect. The crossover number 
A4c of the (T:A) base pairs, where crossover between the strong and weak distance dependence 
occurs, reflects the ratio of intra- and inter-strand neighboring base-base couplings. 

PACS numbers: 87.14.Gg,72.20.Ee,72.25.-b 



In recent years, charge migration in DNA has attracted 
considerable interest among the physics, chemistry, and 
biology communities. Charge transfer in DNA is impor- 
tant for fmictioning of molecular electronic devices 1] as 
well as in understanding the DNA oxidative damage and 
repair Additionally, DNA offers a platform for fun- 
damental physical understanding of systems in the nano- 
scale. It has been a long-standing problem to under- 
stand whether the charge transfer in DNA occurs via a 
unistep coherent superexchange process or a multi-step 
incoherent thermally- induced hopping process 0. In a 
one-strand unistep model the transfer rate is exponential 
and is strongly distance dependent 0, 0, @ . The multi- 
step hopping model on the other hand, predicts a weak 
dependence on the distance. Both of these ideas have re- 
ceived experimental supports 0jl3- Recent experiments 
have shown that the sequence of base-pairs may account 
for the transition between the strong and weak distance 
dependence of charge rates in DNA, but the underlying 
mechanism is not yet clear 0| . In Ref. , the transfer 
rate through a DNA of sequence (G:C)(T:A)a< (G:C)3 
was measured for different A4. The charge transfer 
showed a strong distance dependence when JVl < 3, but 
almost no distance dependence for > 3. To explain 
this distance-dependence crossover at = 3, a combi- 
nation of coherent superexchange and a hopping mech- 
anism (incoherent) - the variable-range hopping model, 
was proposed to allow for a transition between these two 
regimes |^ 0|. In the former process, the donor and 
the acceptor of the charge are coupled to the bridge of 
higher energy, without any chance of intermediate relax- 
ation. The charge remains in a quantum state over the 
bridge that works as a tunneling barrier. In the hopping 
process, relaxation is introduced into each site and the 
charge loses its coherence (phase) when it reaches a site. 
A population parameter for each site was necessary to 
describe the distribution of the charge over the bridges. 

In this letter, we demonstrate that the experimentally 
observed distance-dependence crossover can as well be 
explained by a simple two-strand superexchange model. 
The almost zero distance dependence of the charge trans- 
fer at a long chain is shown to be a result of the inter- 
strand coupling in the DNA. In this model, the system 



still remains coherent and the charge transfer occurs in 
a unistep way. 

We consider a DNA duplex chain of N Watson- 
Crick base pairs connected to four semi-infinite one- 
dimensional (ID) electrodes with one for each end of the 
first and the second strand as illustrated in Fig. ^ The 
tight-binding Hamiltonian of the system is 
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Here cl^ (djj) is the creation operator of holes in the 
first (second) strand on site n of the DNA chain (for 
1 < n < N), the left electrodes (n < 0), and the right 
electrodes {n > N + 1). The on-site energy of site n in 
the first (second) strand is denoted by e„ (wn), which is 
equal to the highest occupied molecular orbit (HOMO) 
energy of the base on this site in the DNA chain and 
the center of conduction band in the electrodes. The 
coupling parameter of the first (second) strand tri,n+i 
{hn^n+i) is equal to the intra-strand coupling parameter 
td between neighboring sites n and n -f 1 of the DNA 
for 1 < n < — 1, one- fourth of the conduction band- 
width in the electrodes t„i for n < — 1 and n > iV -I- 1, 
and the coupling strength tdm between the electrodes and 
the DNA strands for n — and n = N. The inter-strand 
coupling between sites in the same Watson-Crick base 
pair is described by A„. The factor 2 multiplied to each 
sum in Eq. ^ arises from the spin degeneracy. 

In what follows, we have studied the intra-molecular 
hole transfer property along the DNA duplex chain after 
charges are injected (optically or electrically) into the 
base on site 1 of the first strand. To facilitate our cal- 
culation, we connect one virtual electrode to the left end 
of each DNA strand as the injector and another to the 
right end as the drain for holes. To minimize the contact 
effect, we assume a strong coupling (of coupling param- 
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FIG. 1: Schematic illustration of the system. The first strand 
(filled circle) has a DNA base sequence G(T)aiGGG and the 
second strand (empty circle) a sequence C(A)aiCCC. The 
four gray areas indicate the four virtual electrodes connected 
to the DNA chain. Current is injected into the first strand 
through the left electrode Li and measured at the right elec- 
trode Ri. 



eter ioi = tN,N+i = ^oi = hN,N+i = tdm = 1-5 eV) be- 
tween the electrodes and the sites at the ends of the DNA 
strands, and choose a band width (Atm) in the electrodes 



dm 



such that the optimal injection condition td x t„j = t 
[T^ is satisfied. Our result is independent of the choice of 
the value of tdm once it is much larger than the coupling 
parameter between the sites inside the DNA. In this case, 
the added electrodes does not become a bottleneck of the 
system for the charge transfer and the calculated result 
predominantly reflects the properties of the DNA chain. 

The transport properties are evaluated by the transfer 
matrix method |llL Il2l | . For an open system, the secular 
equation is expressed as a group of infinite number of 
equations of the form 



h„^l,n^n-l + {Un - £^)$„ 



A„$„ + 1 



1*. 
1$. 



with 'i'n {^7i) the wave function of the first (second) 
strand on site n. The wave functions of the sites n + 1 
and n are related to those of the sites n and rt — 1 by a 
transfer matrix M, 
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The transmission can be calculated by assuming the 
plane waves propagating in the electrodes. Here we are 
interested in the case where only holes are injected from 
electrode Li to the first strand. The hole wave functions 



in the Li electrode is ^f,, = (Ae'kLua _^ ^^-ikLua s^ < q) 
and in the Ri electrode = Ce''^""" (n > N + 1). The 
distance between two neighboring bases along any DNA 
strand is a = 3.4 A. Using Eq. (Q), we express the output 
wave amplitude C in terms of the input wave amplitude 
A and evaluate the transmission to Ri electrode as 



TiE) 



ICpsin(fcfla) 
|Ap sin(A;La) ' 



We have chosen the normalized incident amplitude to be 
A= 1/V|sin(fcia)|. 

To evaluate the transfer rate or current of a charge 
(hole) from the donor at the left-end site to the acceptor 
at the right-end site of the first strand, we need to know 
the chemical potential at each end. In the experiment 
of Ref. 9], a hole was injected to the left-end site. This 
means that the left chemical potential is approximately 
the on-site energy of this site while the right one is less. 
During the charge transfer process, the hole may retain 
the same energy if no inelastic scattering occurs or lose 
energy via the electron-phonon scattering or other inelas- 
tic collisions. Here we do not deal with these inelastic 
scattering mechanisms explicitly but analyze two limit- 
ing situations, between which the real charge transfer 
process occurs. Since our results for the distance depen- 
dence of the transfer rate from the two limits converge 
(see below), we conclude that our results are reliable. 

In the first limit, we assume that there is no inelastic 
scattering involved and the hole energy is conserved dur- 
ing the transfer process. The transfer rate is proportional 
to the conductance of the system at equilibrium. For a 
small electric potential difference fcsT/e, the current is 



I dET{E)[l-f{E)]f{E). 



(2) 



Here f{E) — 1/ exp[(£' — ^)/kBT] is the Fermi function. 
The room temperature T = 300 K is assumed and the 
on-site energy of site 1 in the first strand is used as the 
chemical potential fi. In the second limit, we assume that 
the hole can lose energy freely during the process, and the 
transfer rate is proportional to the total current via all 
channels of energies below the hole's initial energy. This 
corresponds to an infinitely low chemical potential at the 
right electrode and the current is 



I=—j dET{E)f{E). 



(3) 



We now calculate the distance dependence of the trans- 
fer rate using Eq. (O and Eq. © in a DNA duplex, where 
the first strand has the base sequence G{T) mGGG as in 
the experiment of Ref. . The HOMO energies for bases 
G, C, T, A, are Eg = T-'^^^Ec = 8.87, Et = 9.14, and 
Ea = 8.24 eV respectively [l^. A uniform intra-strand 
hopping parameter tn,n+i = ^n,n+i = td (l < n < N -1) 
and a uniform inter-strand hopping parameter A„ = 
(1 < 71 < N) between any two neighboring bases in the 
DNA are used. 
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First we switch off the inter-strand couphng and cal- 
culate the dependence of the current / on as shown in 
Fig. El (a), for different values of the intra-strand couphng 
parameter td- We find an exponential dependence of the 
current 

/ = /-^ oc e-^^°. (4) 

We then extract the values of f3 for different td and plot 
in Fig. EIc) as /3 vesus \n{td) calculated via Eq. (|3J). The 
curves are almost linear, very similar to the results of 
Eq. (0), and converge to the approximate formula Pi (HIS 



This is the well-known ID superexchange result in the 
literature and has been derived in many different ways. 
This agreement confirms the validity of our model. 

In the next step, we fix td and switch on the inter- 
strand coupling by varying Xd- The result is displayed in 
Fig. |2Ib) where we choose td = 0.5 eV and plot / ver- 
sus for a series of A^. Note that the charge transfer 
occurs via 7r-electrons and generally < td For 
finite Xd, the current drops exponentially with increasing 
Ai for small Ai and then becomes almost flat with oscil- 
lations around a limiting current I°° for large A4. The 
crossover number A4c depends on the strength of the 
inter-strand coupling parameter. The weaker the inter- 
strand coupling is, the bigger the A4c- The dependence 
of /°° on Ad is approximately illustrated in Fig. I2d), 
where the normalized current I^^/I^ of the DNA chain 
at = 10 is plotted versus In(Ad). Again, two almost 
identical straight lines are found corresponding to the two 
limiting situations based on Eq. ^ and Eq. Q and can 
be approximately expressed as 

ln(/i7/i) = 5.7 + 3.91n(Ad). (6) 

From Eqs. Q-©, we estimate the ratio of inter- and 
intra-strand coupling from the crossover number Mc- 
Since the environment can change Xd/td, we predict that 
the transition number may vary and be different from 3 
when the experimental environment changes. 

Calculating the current / before and after adding a 
(T:A) base pair at site n with zero or nonzero inter- 
strand coupling A„, we find that the distance-dependence 
crossover has a topological origin, viz., from the ID chain 
charge transport to a partly two-dimensional (2D) net- 
work. When a new (T:A) base pair is inserted into 
the DNA chain, a new superexchange channel is opened 
through its inter-strand coupling and the correspond- 
ing contribution exactly compensates the loss of charge 
transfer rate that would incur because of an extra barrier 
to the existing channels. 

In Fig. 121 we fit the Ai dependence of the charge trans- 
fer rate observed in Ref. using intra- and inter-strand 
couphng parameters td = 0.52 eV and A^ — 0.07 eV re- 
spectively. Eq. ||2J) is employed in the calculation. The 




FIG. 2: (a) Current / versus M for ta — 
0.1,0.2,0.3,0.4,0.5,0.6 eV (from lower to upper curves) 
for zero inter-strand coupling. The displayed results are from 
Eq. @ and identical results are obtained from Eq. @ in 
all the panels, (b) Same as in (a) at fixed td = 0.5 eV but 
for Xd = 0, 5, 20, 40, 80, 100 meV corresponding to curves 
counted from the bottom, (c) The f3 value calculated from 
the slope of the lines in (a) versus Intd. (d) ln(/^''//^), where 
J"^ is the current for a chain with M (A:T) base pairs, 
versus InA^j. The unit of td and is eV. 



agreement between the experimental and theoretical re- 
sults are very good except that a small oscillation is vis- 
ible in the theoretical result near A4c- This oscillation 
results in the deviation of the empty circle from the filled 
circle at = 4. When Eq. Q is used, similar result is 
obtained but with a stronger oscillation. The oscillations 
reflect the fact that we have treated the system as a co- 
herent system by neglecting the dephasing effect from the 
environment and the relaxation process from phonons. 

To get a clear picture of the process, we plot as inset 
in Fig. 3, the transmission T as a function of the hole en- 
ergy E for systems with A^ = 1, 2, 3, and 7 in an energy 
range near and below the G base HOMO energy Eg- In 
the T spectrum, each peak represents a transport channel 
and there are more fine structures or peaks when more 
base pairs are added to the system. When A4 varies from 
1 to 3, the ID chain transport dominates and only one 
principal transmission peak is important. The principal 
peak shifts when Ai varies due to the shift of energy of 
the channel; its height drops rapidly leading to an expo- 
nential decrease of charge transfer rate. If we add more 
(T:A) base pairs to the DNA duplex, the principal T peak 
drops to a level comparable to that of other peaks and 
results in a corssover from ID chain transport to 2D net- 
work transport. In the absence of any inelastic scattering 
the charge transfer rate versus A4 oscillates as a result of 
the energy shift of the transport channels and the energy 
conservation of the charge. With the assistance of the 
phonon, however, the charge can use channels of energy 
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FIG. 3: Normalized transfer rate measured in Ref. j^] (filled 
circle) and theoretical fit using this model (open circle), 
log(/'^//'^), are plotted as functions of M (T:A) base pairs 
between the (G:C) and the triple (G:C) base pairs. Inset: 
The corresponding transmission T versus energy E for A4 = 1 
(solid line), 2 (dotted line), 3 (dot-dashed line), and 7 (dashed 
line) . 



different from its initial energy and phonons may play an 
important role in assisting the charge transfer. 

Our intra-strand coupling parameter td = 0.52 eV used 
to fit the measurement is consistant with the GG cou- 
pling parameter extracted from a direct I-V measurement 
through a DNA of 30 (G:C) base pairs This fit pa- 
rameter is much larger than the ab initio values [13, uM ■ 
The reason of this disagreement is not yet clear but it 



may be related to the electron-phonon interaction. This 
interaction may also affect the temperature dependence 
of the charge transfer through the system In this 

work, we do not treat the electron-phonon interaction in 
detail but we expect that these interaction can determine 
the position of the real current between our two current 
limits [Eq. (|2J) and Eq. ©J discussed above. The de- 
phasing effect due to the environment can also result in 
a weak distance dependence for a ID long bridge system 
l^i, which is a different mechanism from what we pro- 
pose here. The dephasing effect exists in a real system 
and can help damp the oscillation of the current observed 
in Figs. 13 and El 

In summary, we propose a new mechanism for the 
charge transfer through a DNA duplex chain. It is dif- 
ferent from the previously proposed thermally-induced 
hopping mechanism in explaining the observed weak dis- 
tance dependence when the number of (T:A) base pairs 
between (G:C) base pairs is larger than A^c, in that we 
treat the system fully quantum mechanically and em- 
phasize the importance of inter-strand coupling between 
the two strands of the DNA duplex. We found that the 
series of (T:A) base pairs in long (G:C)(T:A)>i (G:C)3 
DNA duplex chains is still a quantum tunneling barrier. 
The holes in the left (G:C) base pair tunnel through this 
two-strand network superexchange barrier instead of one- 
strand chain superexchange barrier, to the right triple 
(G:C) base pairs. 
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